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Abstract. In recent years higher-dimensional black holes have attracted much interest because 
of various developments in gravity and high energy physics. But whereas higher-dimensional 
charged static (Tangherlini) and uncharged rotating (Myers-Perry) black holes were found long 
ago, black hole solutions of Einstein-Maxwell theory, are not yet known in closed form in more 
than 4 dimensions, when both electric charge and rotation are present. Here we therefore 
study these solutions and those of Einstein-Maxwell-dilaton theory, by using numerical and 
perturbative methods, and by exploiting the existence of spacetime symmetries. The properties 
of these black holes reveal new interesting features, not seen in D = 4. For instance, unlike the 
D — 4 Kerr-Newman solution, they possess a non-constant gyromagnetic factor. 



1. Introduction 

In Einstein-Maxwell (EM) theory in 3 -|- 1 dimensions the unique family of stationary 
asymptotically flat black holes with non-degenerate horizon comprises the rotating Kerr-Newman 
(KN) and Kerr black holes and the static Reissner-Nordstrom (RN) and Schwarzschild black 
holes. EM black holes are uniquely determined by their mass, their electric and magnetic charge, 
and their angular momentum (see [T] and references therein). 

The generalization of black hole solutions to higher dimensions was pioneered by Tangherlini 
[2] for static vacuum black holes, and by Myers and Perry (MP) [3] for stationary vacuum 
black holes. Whereas Myers and Perry [3J also obtained charged static black holes in higher 
dimensional EM theory, higher dimensional charged rotating black holes have not yet been 
obtained in closed form in pure EM theory, although such black holes are known for some low 
energy effective actions related to string theory (including additional fields, though) [H [5l [6l [7] . 

Higher dimensional black holes received much interest in recent years also with the advent of 
brane-world theories, raising the possibility of direct observation in future high energy colliders 
[8], which makes the theoretical study of these higher dimensional objects particularly appealing. 

However, finding black holes solutions in higher dimensions is a difficult task due to the size 
and complexity of the equations. Moreover, contrary to what happens in D = 4 the topology of 
the horizon is not unique in higher dimensions (for instance, ringlike configurations are allowed 
[9]), which makes the problem even harder. For these reasons the number of solutions known 
analytically in closed form is very limited. This forces us to use alternative techniques to study 
these black holes. Here we will briefly analyze such black holes by two methods: the numerical 
method and the perturbative method. 



2. Einstein-Maxwell-Dilaton theory 

We will concentrate on Einstein-Maxwell-Dilaton (EMD) theory in D dimensions, as an example 
how these methods may be applied. The classical EMD action reads 
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where R is the scalar curvature, $ is the dilaton field, F^j^ = dfj,Ai, — d^A^ is the electromagnetic 
field tensor, and is the electromagnetic potential, h is the dilaton coupling constant. 

The equations of motion can be obtained by varying the action with respect to the metric 
g^y^ the dilaton field $ and the gauge potential A^, yielding 
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For /i = the theory reduces to pure EM theory. 

These equations admit an analytical solution when the dilaton coupling constant h takes the 
Kaluza-Klein (KK) value Hrk ffOl E] 
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In order to generate charged rotating solutions for generic values of h we will assume that the 
topology of the horizon is spherical. We will further assume that the black holes are stationary 
and axisymmetric. Stationary axisymmetric black holes in D dimensions possess = \{D—l)/2\ 
independent angular momenta Jj associated with N orthogonal planes of rotation [3|. {[{D—1) /2] 
denotes the integer part of (D — 1)/2, corresponding to the rank of the rotation group SO{D — l).) 
These black hole solutions then fall into two classes: even-Z) and odd-D-solutions. 

When all N angular momenta have equal-magnitude, the symmetry of these black holes is 
strongly enhanced. In fact, in odd dimensions, the symmetry then increases from R x U{i)^ 
to R X U{N), thus changing the problem from cohomogeneity-A^ to cohomogeneity-1. As a 
consequence the original system of partial differential equations in variables then reduces to 
a system of ordinary differential equations (ODE's). 

3. Numerical and perturbative methods 

In what follows we will focus on the case where D is odd, all angular momenta have the same 
magnitude, and the topology of the horizon is that of a sphere. That case is tractable numerically 
and perturbatively. Under those assumptions the metric, the gauge potential and the dilaton 
field can be parametrized with their angular dependence explicitly known 
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where = 0, 9i € [0,7r/2] for i = 1,...,A^-1, 6^ = vr/2, ipk G [0, 27r] for = 1, . . . , A^, 
and £k = ±1 denotes the sense of rotation in the k-th orthogonal plane of rotation. The 
spacetime possesses + 1 commuting Killing vectors, ^ = and ry(fc) = 5^^., for A; = 1, . . . , A^. 
This parametrization employs isotropic coordinates and is appropriate for numerical purposes 
\12\ I13j. However, for perturbative purposes Boyer-Lindquist coordinates simplify the algebra 
[a US]. 

In order that the solutions correspond to black holes they must possess an event horizon 
located at r = rn, characterized by the condition /(rn) = [TE]. Such an event horizon is in 
fact a Killing horizon since the Killing vector 
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is null at and orthorgonal to the horizon, Q representing the constant horizon angular velocity. 

The system of ODE's resulting when Eqs. ([5ll6|) are substituted in Eqs. (I2]l3]) must be 
supplemented with the appropriate set of boundary conditions. This is obtained by imposing 
that the horizon is regular and the spacetime is asymptotically flat [111 \T2\ I13j . Then the 
system of ODE's is solved by employing a numerical solver (e.g., COLSYS), providing accurate 
results. The numerical parameters are {/i, rn, f^, Q}, with Q being the electric charge. Varying 
these parameters we obtain families of solutions, with their corresponding values of the physical 
properties like the mass M, the angular momentum J, etc. An important information the 
numerical method provides is the domain of existence of black holes in higher dimensions. That 
domain consists of the region in the parameter space bounded by the extremal solutions. In 
Figure 1 (left) we exhibit the domain of existence of EMD black holes in terms of the scaled 
electric charge Q/M and angular momentum J /M^^~'^^/^^~^\ Clearly the domain depends 
both on the dilaton coupling constant h and the dimension D. 
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Figure 1. Left: Domain of existence of EMD black holes with equal-magnitude angular 
momenta in D = 5, 7 and 9 dimensions. Right: Comparison of the gyromagnetic factor g of 
D = 5 extremal EMD black holes obtained with numerical and perturbative (cyan) methods. 

The accuracy of the numerical method may be tested by means of the analytical KK solutions, 
as well as by exact relations (like the Smarr mass formula) these black holes satisfy [15j . 

As a complementary method one may employ perturbations to generate black hole solutions 
in higher dimensions [TH [15l [l7] . The advantage of this method is that it produces analytical 
formulae, which the numerical method does not. However, the range of validity of these formulae 
is limited by the accuracy of the order of the perturbations. 



A good perturbative parameter to produce charged rotating black hole solutions in higher 
dimensions is the electric charge. One starts from the MP solution [3j and expands the metric, 
the gauge potential and the dilaton in series expansions in the electric charge around the MP 
solution (see [I^ for details). Those series are substituted in Eqs. dSJlS]) and solved order by 
order in the electric charge. For lower orders the equations are easy to solve but as the order of 
the perturbations increases the equations become more and more difficult to solve. 

An important point when applying the perturbative method is the way the integration 
constants are fixed. One has to impose regularity of the solutions at the horizon together 
with the condition of asymptotic flatness. The Smarr mass formula is a useful means to check 
the accuracy of the solutions. 

When the perturbative method is applied to extremal solutions the formulae simplify very 
much. For instance, one obtains the following formula of the gyromagnetic factor g to second 
order 

showing that the gyromagnetic factor is not constant for pure EM solutions (/i = 0). 
4. Comparison of the two methods 

Although both methods have a different range of applicability there is a region of the parameter 
space where both are valid, namely, when the electric charge is small. In that region one 
can check the consistency of the methods and their level of agreement [15] ■ This is presented in 
Figure 1 (right) where the perturbative formula Eq. ([5]) for extremal solutions is compared to the 
corresponding non-perturbative numerical results for D = 5. We observe that for \Q/M\ < 0.2 
both approches are in good agreement. In fact, the perturbative method may be modified to 
analyze also slowly rotating solutions, using as the unperturbated initial solution the Tangherlini 
solution and using the angular momentum as the perturbative parameter. 

To summarize, these two methods are complementary and provide us with a powerful tool 
to analyze black holes in higher dimensions when analytical solutions are not available. 
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